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Abstract. We propose a model for non isothermal ferromagnetic phase transition based 
on a phase field approach, in which the phase parameter is related but not identified with 
the magnetization. The magnetization is split in a paramagnetic and in a ferromagnetic 
contribution, dependent on a scalar phase parameter and identically null above the Curie 
temperature. The dynamics of the magnetization below the Curie temperature is gov- 
erned by the order parameter evolution equation and by a Landau-Lifshitz type equation 
for the magnetization vector. In the simple situation of a uniaxial magnet it is shown 
how the order parameter dynamics reproduces the hysteresis effect of the magnetization. 

1. Introduction 

The peculiar feature of ferromagnetic materials is the behaviour of the magnetiza- 
tion vector below a characteristic value 6c of temperature, named the Curie temperature. 
Firstly, at temperatures 6 < 9c, a, non zero value Mo(^) of the magnetization is stable even 
at zero external magnetic fields; this magnetization is said the spontaneous magnetization. 
On the contrary, in the paramagnetic regime, that is for 9 > 9c, the magnetization vanishes 
at zero external field. Moreover, the way in which external fields influence the magnetiza- 
tion vector differs in the two cases. While in the paramagnetic state the magnetization at 
a point is a function of the magnetic fields at that point (with a proportionality relation 
at sufficiently low fields) , in the ferromagnetic state the magnetization shows well known 
hysteresis phenomena [1]. In the ferromagnetic materials the external fields influence 
primarily the direc tion of the magnetization vector, and a model for the magnetization 
evolution in saturation conditions has been proposed long ago by Landau and Lifshitz 
[2, 3]. 

These qualitative differences in the behaviour of the magnetization on the two sides of 
the Curie temperature can be understood in the framework of the phase transitions. The 
clarification of this issue is a fundamental contribution of the Landau theory of phase 
transitions [4, 5]. The approach of Landau is based on the concept of an order parameter, 
which is a physical (macroscopical) observable quantity whose behaviour is able to account 
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for the microscopical change of structure which generally characterizes phase transitions. 
Strictly speaking, the Landau theory accounts for second order (or continuous) phase 
transitions, which, according to Landau, are properly understood in terms of symmetry 
breaking. So, in the case of ferromagnetic transition, for example, the order parameter of 
the Landau theory is the magnetization vector: the transition manifests itself as rotational 
symmetry breaking due to the set up of a non-zero magnetization 

Mo under otherwise isotropic conditions, that is in zero external field. Perhaps, it is not 
superfiuous to point out that the value of the order parameter is not in itself an indicator 
of the phase of the material, except in the particular case of null external field. 
Along these lines, in [6] a three dimensional evolutive model is proposed, using the whole 
magnetization vector with the order parameter, ruled by a vectorial Ginzburg-Landau 
time-dependent equation. 

In this paper we propose a model different from the original Landau setting and closer to a 
general phase field approach. That is, we introduce a scalar phase field, which (unlike the 
magnetization vector) vanishes above the Curie temperature, even in presence of external 
fields. The order parameter vanishes in a continuous way when the Curie temperature 
is approached, according to the second order character of the transition. The magneti- 
zation will be decomposed in two contributions, one which is of paramagnetic character, 
namely it is a direct function of the external field, and an other one, depending also on 
the order parameter, which is considered as an independent field with its own evolution 
equation (the time-dependent Ginzburg-Landau equation). The sense of this separation 
is not, of course, that of a physical distinction between two different sources of magneti- 
zation, as well as, for example, the two fiuid theory of superfiuidity [7] is not the theory 
of a mixture of fiuids. The two magnetization contributions refiect the different way in 
which the magnetization evolves under the external field in the ferromagnetic and in the 
paramagnetic regime. In particular, the hysteresis phenomena manifesti ng themselves 
in the ferromagnetic phase are determined by the coupling with the phase field evolution 
equation. Under this respect, the model draws on the internal variable models, in which 
a history-dependent constitutive equation for a physical quantity (in our case the mag- 
netization) is obtained through the coupling with an internal variable obeying its own 
differential evolution equation ([8]). Nevertheless, in this model, the phase field is not a 
mere internal variable, nor its evolution equation is a constitutive one. We assume that 
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the phase field equation is a true balance equation associated with its own power balance 
[9]. 

The phase field of this model is related to the modulus of the spontaneous magnetization 
(nevertheless it is infiuenced by the external field); loosely speaking, we can say that its 
physical meaning is related to the microscopic order set up by the microscopic exchange 
interactions, responsible of spontaneous magnetization. 

The model we are proposing couples the Maxwell equations for the electromagnetic field 
with a scalar time-dependent Ginzburg-Landau equation for the order parameter and a 
heat balance equation for the temperature. Two constitutive equations define the rela- 
tion between the magnetization and the order parameter: the magnetization is split in 
a paramagnetic and a ferromagnetic contribution and the direction of the last one, in 
the tridimensional case, is ruled by a Landau-Lifshitz-Gilbert equation [2, 10], suitably 
modified in order to describe the non saturated regime. The constitutive choices of the 
model are verified to be consistent with the second law of thermodynamics in the form of 
the Clausius-Duhem inequality. 

Finally, we present a particular case of the model when the direction of the magnetization 
is fixed, as in the case of uniaxial ferromagnets (one-dimensional model). We show in this 
situation how the dynamics of the order parameter in the ferromagnetic phase gives the 
usual shape of the hysteresis cycle of the magnetization. 

2. Three dimensional model 

Let us consider a ferromagnetic material occupying a bounded domain 1] C M^. Denot- 
ing by E, H, D, B the electric field, the magnetic field, the electric displacement and the 
magnetic induction, the behavior of the material is ruled by Maxwell's equations 

(1) VxE = -B, VxH = D + J, 

(2) V-B = 0, V-D = p, 

where J is the current density and p is the charge density. We assume the constitutive 
equations 

(3) D = £E, B = hqR + M, J = crE, 

where e,fj,Q,a are respectively the dielectric constant, the magnetic permeability and the 
conductivity, while M denotes the magnetization. As known, in paramagnetic materials. 
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the magnetization is a function of the magnetic field H. On the contrary, ferromagnetic 
systems are characterized by a time-non local relation between magnetization and the 
magnetic field. Therefore, in order to describe paramagnetic-ferromagnetic transitions, 
we write 

(4) B = /xH + M, /x = /x(H,e). 

This amounts to split the magnetization as M = (// — //o)H + M, where M is the part of 
the magnetization whose value at a given time cannot be expressed as function of the field 
H at the same time. This contribution exists only in the ferromagnetic phase and it is 
history-dependent, so M = in the paramagnetic state and M 7^ in the ferromagnetic 
regime. The model we propose is set in the general context of the Ginzburg-Landau 
theory by defining, as order parameter, a scalar phase variable </? such that (p — in 
the paramagnetic phase and </? > in the ferromagnetic state. As a consequence the 
magnetization is related to (p. More precisely, we assume that 

(5) M = M{ip, e)m, M{ip, e)>0 

where m is a unit versor and the modulus M{ip, 6) depends also on the temperature. 

The evolution of the phase ip is given by the Ginzburg-Landau equation typical of 
phase transition models. We introduce the classical potentials describing second-order 
phase transitions ([9]) 

and assume the following equation 

(6) Tip = ^Aip - ecF'iip) - eG'iif) - A{ip, e)i{ec - ^)h • m 

where T, arc positive constants, A is a generic function whose definition will be specified 
later, and / is the unit step function, i.e. 

if x < 

1 if X > 



I{x) 



Notice that for large values of the temperature, 9 > 6c, last term of (6) vanishes and the 
function 

w{ip) = ecF{ip) + eG{ip) 

admits the minimum value (/? = which characterizes the paramagnetic phase. 
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The evolution equation for ip should preserve the defining condition 97 > , which is not 
automatic in eq. (6). So this is a further constraint which has to be enforced, for example, 
in any numerical solution of the equation, and amounts to add a singular contribution in 
the potential function F{(p) such that F{ip) = +00 for ip < 0. 

Concerning the evolution of the direction of the magnetization, we assume that the 
unit versor m satisfies the Landau-Lifshitz equation ( [2] ) 

(7) (p^m = -7(/7^m x H - Am x (m x H), 7, A > 0. 

Notice that 

(/^'|m| = (7V^ + A')^|mxH|. 

As a consequence, when ip approaches zero, the direction of m moves toward the direction 
of H. Moreover, by multiplying (7) by m we obtain 

(/7^m • m = 0, 

which is consistent with the condition |m(x, t) \ — 1 for any t > 0, provided that |m(x, 0) | = 
1. 

Like other models of phase transitions (see [9]), equation (6) can be interpreted as a 
balance law of the order structure. Indeed it can be written in the form 



A; = V-p, 



where 



k = rip + ecF'{ip) + eG'{ip) + A{ip,e)i{ec-e)H-ni, 

p = ^Vv?. 

This formulation allows us to define the internal power related to the phase variable as 



= + p • V0 



1 



Tip' + OcF^ip) + eG{ip) + A{ip, e)i{ec - e)ipin ■ h + — v<^ ■ 



From (4) and (5) we deduce the relation 



B = /^H + ( t ^ + ^ • h1 H + f ^0 + ) m + Miip, e)rn. 



de dU J Kdip"^ 09 
Hence the electromagnetic power 

T'ei = B H + D E + crE 
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can be written as 

+£E • E + (tE^ 

We denote by e the internal energy and h the thermal power. The first law of thermody- 
namics reads 

(8) e^Vei + V^ + h. 
where h satisfies the thermal balance law 

(9) /i = - V • q + r, 

and q, r are respectively the heat fiux and the heat source. 

In order to prove the consistence of the model with the second law of thermodynamics 
we look for the constitutive relations for the entropy function rj and the heat fiux q that 
ensure the fulfillment of Clausius-Duhem inequality 

..-v.(|).:. 



Thermal balance law (9) yields 



9ri>^-V9 + h. 
u 



i^ + vO<Vei + v^-^- ve. 



Hence, by introducing the free energy ip — e — Orj, the previous inequality leads to 

q 

e 

By substituting the expressions of the powers, we obtain 

i. + ,e < ;.H.H+(|^ + |^.H)H^+(^0 + ^^)m.H + M(^,^)m.H 

+£E • E + (te^ + Tip^ + ecF(ip) + eG(ip) + A(ip, e)i{ec - e)ipxa. ■ h 

By means of (7), we deduce 

+AM((^, e)Lp-^\m. X H|2 + £E • E + aE^ + r^^ + 9cF{lp) + 9G{lp) 

e)i{ec - e)(pm ■ h + \\/(fi • v^? - § • ve. 
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We assume that the free energy depends on the variables (99, 6*, E, H), so that 
the previous inequahty yields 

de do' 



'di) 1 




+ 


" dip 




+ 


dip 


dH. ' 



Ya.-ll-e^F'{ip)-eG'{ip) 
dM 



+ 77 - - • H 



e + 



dip 



9E 



-£E 



E 



H 



(10) < XM{ip, e)ip-'\m. X H|^ + aW + Tip^-^- V9. 

The previous inequality is fulfilled if we choose the constitutive relations 

dM 

q^-ko(e)ve ko[e)>o. 

Usual arguments of thermodynamics based on the arbitrariness of {(p, V(p, 6, E, H) lead 
to the following expressions of the free energy and entropy 



^ 89 89 89 



m H 



Substitution into (10) yields 



dM 
~d9 



m H. 



\M{ip, e)^-'\m X H|' + (tE^ + Tip' + 



HO) 



\V0\^ > 0, 



which guarantees that Clausius-Duhem inequality is satisfied. 

The evolution equation for the temperature follows from the thermal balance law (9), 
by substituting the expression of h deduced from the first law (8). Since the internal 
energy is written as 



^-^2 , ..tt2 ^ 



where 



e = ^ + 0rj^eoi9) + ^-E' + i^li' + -^\Wip\'' 

1^ - 9^\ H • dH + 9,F{^) + 9^m ■ H, 



eo(9)^iPo(0)-0m, 
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substitution into (8) yields 

(11) -Tif'^ - aE'^ - \M{ip,0)ip-^\mxU\'^. 
Hence the temperature satisfies the equation 

(12) h^-V ■[ko{e)ve]+r. 

In this model we assume the following constitutive equation 

A{<p, 9)^9- 9c. 

This choice leads to a continuous temperature dependence for the modulus of M, namely 

M^^{9- 9c)-, 

where the subscript — denotes the negative part of a function, i.e. /_ = max{— /, 0}. 
Therefore the Ginzburg-Landau equation for the phase field reads 

Tip = \Aip - 9cF\ip) - 9G\ip) + {9- ^e)-H • m. 

3. One dimensional model 

In this section we will consider a one-dimcnsional model, obtained by assuming that 
the magnetic and electric fields have constant and orthogonal directions, say y, z, and 
that the components of the unknown fields depend only by the variable x, namely 

E = E{x)\i, H = H{x)l 

In the description of uniaxial ferromagncts, wc modify the definition of the order param- 
eter, by requiring that (/? 7^ in the ferromagnetic phase and (/? = in the paramagnetic 
state. Therefore (p is allowed to take negative values and the vector m is defined as 

m = sign{(p)j. 

We assume the constitutive equations 

(13) M{ip,9) = \^\{9-9c)-. 
Therefore 

M^cp{9-9c)-j 
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and equations (l)-(4) imply 

eE = d^H-aE 

^H+^H9+^HH = d,E-<p{9-9c)_ + <pI{9c-9)e. 

The evolution of the magnetization (p is governed by the Ginzburg-Landau equation 

1 

Finally, the evolution equation for the temperature is deduced by (11) and (12) 



(14) Tip = -d,,^ - 9cF'{ip) - 9G'{ip) + {9- 9c)-H. 



_ 0G^^) + g^HH - 9j^ [I{9c - 9)^H] 



e'oi0) + 9j^ll.dil 
(15) = Tip'^ + aE^ - d^q + r. 



A constitutive equation of ii{9, H) has to be given. For example, in the classical Landau 
model of ferromagnetism, the (total) magnetization M as a function of the magnetic field 
and the temperature is given by 

(16) hoM^ + ao{9 - 9c)M -H^O, 
from which it is obtained the permeability aX 9 > 9c 

(17) ii-1 



/io hoM{H,9f + ao{9-9c)' 

For H — this equation provides the well known Curie- Weiss law for the susceptibility, 

Xo oc 9 > 9c. 

For H ^ the resultant permeability is a regular function of the temperature. In our 
model, it is required that, whatever the constitutive relation for /i is taken, the integral 



exists finite. We observe that for if — > oo, d'^ii/d9'^ is expected to tend at zero for satu- 
ration reasons, while, in this respect, the permeability resulting from (17) is reasonable 
only at small fields. Moreover, the function eo(^) has to satisfy eo(^) -|- 9J{9, H) > for 
every 9 and H to have a standard parabolic heat equation. 

We see that equation (14) is able to recover the hysteresis diagram typical of the phe- 
nomenon of ferromagnetism. To this purpose, we will consider a spatially homogeneous 
material in isothermal conditions, with 9 <9c. Then (14) reduces to 

(18) Tip = -9cF'{ip) - 9G'{ip) -{9- 9c)H. 
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Moreover, from (4) and (5) we obtain 

(19) B^iiH -{e- ec)ip. 

Here we assume /i as approximately if-independent in the considered range of the mag- 
netic field. If H is a known function of time, equations (18)-(19) allow us to obtain 
the B — H diagram. In particular ii H = HQsm{ujt), t G [0,T] and the initial condi- 
tion is (/9(0) = 0, we deduce the following hysteresis diagrams for different values of the 
temperature. 




Figure 1. B — H diagram with the numerical constants 9 — 0.9, 9c — 
1, a; = 7r, r = 0.01, Hq ^ A, = 0.01, T = 2.5. 




Figure 2. B — H diagram with the numerical constants 9 — 0.5, 9c — 
1, co^n, T = 0.01, i/o = 4, = 0.01, 7 = 2.5. 
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Figure 3. B — H diagram with the numerical constants 9 — 0.1, 9c — 
1, uj^n, T = 0.01, Ho ^4:, n = Q.Ql, T = 2.5. 
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